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ABSTRACT 

Context. The importance of the physical environment in the evolution of newly formed low-mass stars remains an open question. In 
particular, radiation from nearby more massive stars may affect both the physical and chemical structure of these kinds of young stars. 
Aims. To constrain the physical characteristics of a group of embedded low-mass protostars in Corona Australis in the vicinity of the 
young luminous Herbig Be star R CrA. 

Methods. Millimetre wavelength maps of molecular line and continuum emission towards the low-mass star forming region IRS7 
near R CrA from the Submillimeter Array (SMA) and Atacama Pathfinder Experiment (APEX) are presented. The maps show the 
distribution of 18 fines from 7 species (H2CO, CH3OH, HC3N, C-C3H2, HCN, CN and SiO) on scales from 3" to 60" (400-8 000 AU). 
Using a set of HtCO lines, we estimate tfie temperatures and column densities in tfie region using botfi LTE and non-LTE metfiods. 
The results are compared with 1-D radiative transfer modelling of the protostellar cores. Tfiese models constrain which properties of 
the central source, protostellar envelope, and surrounding radiation field can give rise to the observed line and continuum emission. 
Results. Most of the H2CO emission from the regions emerges from two elongated (~ 6000 AU long) narrow (< 1 500 AU) ridges 
dominating the emission picked up in both interferometric and single-dish measurements. The temperatures inferred from the H2CO 
lines are no less than ~ 30 K and more likely 50-60 K, and the line emission peaks are olfset by ~ 2 500 AU from the location of 
the embedded protostars. These temperatures can not be explained by the heating from the young stellar objects (YSOs) themselves. 
Irradiation by the nearby Herbig Be star R CrA could, however, explain these high temperatures. The elevated temperatures can in 
turn impact the physical and chemical characteristics of protostars, in particular, lead to enhanced abundances of typical tracers of 
photon dominated regions (PDRs) such as seen in single-dish line surveys of embedded protostars in the region. 

Key words, stars; formation - ISM: individual objects (R CrA) - ISM: molecules - radiative transfer - astrochemistry 



1. Introduction 

Low-mass stars form when dense molecular cloud cores un- 
dergo gravitational collapse. In the very early stages (Class 0; 
|Andre et al.||1993 [2000 j l, the young star is obscured by its en- 
velope of gas and dust. This envelope eventually dissipates, and 
leaves a pre-main-sequence star with a planet-forming disc. One 
of the interesting questions concerning the star formation pro- 
cess is how impoitant the ambient environment is for the out- 
come, e.g. in terms of the physical and chemical properties of 
the emerging protostars and their discs. It is for example well- 
known that high-mass stars in the environment of young stellar 
objects (YSOs) can aff'ect both the chemical and physical evolu- 
tion of the protostars (see e.g. Tielens & Hollenbach 1985a bjl. 
The most well-known examples of young objects aff'ected by 
highly luminous stars are the proplyds in the Orion Nebula, 
which were discovered in Hubble Space Telescope observations 
( O'Dell et al.||19931 ). In these sources, the luminous emission 
from the high-mass stars results in photoevaporation of the pro- 
toplanetary discs ([ Johnstone et al.|1998) ). In this paper, we focus 
on the contribution of external irradiation onto more deeply em- 
bedded sources in the Corona Australis (CrA) region. 



The CrA dark cloud is one of the most nearby active star 
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of the combined SMA 
electronic form at the CDS via 



* Full channel maps (Figs, 
and APEX data are available m 
anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http : //cdsweb . u-strasbg . fr/cgi-bin/qcat? J/A+A/ 



forming regions at a distance of 130 pc (Neuhauser & Forbrich 



2008|). Throug h multi-wavelength studies of the CrA region, 
Peterson et al.| ( |20TT] l found 116 YSO candidates, the vast ma- 



jority being Class II and Class III sources. The region around the 
Herbig Be star R CrA - which is also referred to as the Coronet 
cluster - contains the highest concentration of newly formed and 
young stars (Class I and earlier; Peterson et aI.|20lT)l. R CrA it- 
self is a Herbig Be star (spectral class B 5-B8;|Gray et al.|2006t 
BiboetaI.|1992| with a mass of 3.0 M© (Bibo et al.|1992| l and a 



luminosity between 100 and 166 Lq (the lower value is a lower 
limit from SED estimates with incomplete data in the UV range; 
|Bibo et aL]|1992] l. Several dense molecular cloud cores with 
masses between 2 and 50 Mq were fou nd near the Herbig Be star 



R CrA in C"*0 SEST observations by|Harju et al.|(!l993). The 



IRS7 region is located with its centre between the Herbig Be star 
R CrA and the T Tauri star T CrA ( [Taylor & Storey] 1984| l, -30" 
(a: 4 000 AU) from each of these stars, and harbours a handful of 
Class OA YSOs. 

A summary of continuum point sources in the region around 
R CrA reported in the literature is given in Table [Tl The lo- 



cations of the sources are shown on Spitzer 4.5 \im ( [Peterson 
et aI.||20l"T| l image data, with Herschel 70 (tm (Lindberg et al., 
in prep.) and SMA 1.3 mm (|Chen & Arce||2010|l continuum 



data overplotted in Fig. [T] Using NIR photometry, Taylor & 
|Storey| ( |1984| l found that IRS7 is the most reddened source in 
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Fig. 1. SMA 1.3 mm continuum (lOcr solid contours, cr - 4.7 mJy beam '), Herschel PACS 70 [im continuum (10 Jy dashed 



contours; see also Lindberg et al., in prep.), and Spitzer 4.5 ^m image (greyscale; see also |Peterson et al. [201 Ij . The Herschel 
PACS beam size is 9", and the data consist of two overlapping 5x5 spaxel arrays, with a 9.4"x9.4"spaxel size. The SMA beam 
size (6.0" x2. 9") is shown in the bottom right corner. The plus signs indicate the positions of the YSOs ( Peterson et al.|20lT ) and 
SMM 1 A ( Nutter et al. 2005 ), and the star symbols indicate the positions of R CrA and T CrA ( ,Peterson et al. 201 1 1. These symbols 
will be used throughout the paper. 



the R CrA region, with a visual extinction of more than 25 mag. 
Through VLA 6 cm observations, |Brown| ( 1987) 1 discovered 



that 1RS7 is of binary nature, and separated it into IRS7A and 
IRS7B, two apparently compact sources with a separation of 10- 
15" (2 000 AU). |Brownl(|1987|l also detected a faint continuum 



source north of IRS7A. [Miettinen et al.| ( |2008j ) refeiTed to this 
as Brown 9, and used 3 cm and 6 cm continuum observations 
to produce high-resolution high-sensitivity maps of the region. 
Nutter et al.| (200 5 ) studied the region in submillimetre, and 
detected three peaks within the region around 1RS7: SMM IB 
coinciding with IRS7B, SMM IC coinciding with Brown 9, 
and SMM 1 A south of 1RS7B. Studies of millimetre continuum 
emission showed three compact sources, IRS7B (SMM IB), 

[Chen & Arce|20T0l 



SMM IC, and SMM 2 (Groppi et al.|2007 



Peterson et aLl|2011| l. [Groppi et al.| ^2007| l used Spitzer IRAC 



data, SCUBA 450 ^im and 850 ^im, and SMA 1.1 mm data to 
study the SEDs of 1RS7B (SMAl) and SMM IC (SMA2; not de- 
tected in the Spitzer data). They suggested that 1RS7B is a tran- 
sitional Class 0/1 object, whereas SMM IC likely is a Class 
source. IRS7A shows strong mid-lR emission (Spitzer), but is 



not detected in the (sub-)millimetre. Based on their IR and sub- 
millimetre SEDs, [Peterson et al.| ( [201 1[ ) suggested that 1RS7A 
and SMM 2 are Class I young stellar objects. SMM lA is sug- 
gested to be a prestellar core ( Chen & Arce|2010.) . 



Compared to other typical embedded protostars 1RS7B 
has a fairly rich line spectrum and is of interest as a pos- 
sible astrochemical laboratory. Detailed studies of the chem- 
istry of 1RS7B have been made possible with a new genera- 
tion of southern hemisphere submillimetre telescopes such as 
Atacama Pathfinder Experiment (APEX) ( [Schoier et al. 2006[ l 
and Atacama Submillimeter Telescope Experiment (ASTE) 
(Watanabe et al. 2012i. [Wa tanabe et al.' performed a full line 
scan of IRS7B in the 345 GHz submillimetre window using 
ASTE, and found H2CO and CH3OH, two molecules typical 
for a chemistry dominated by evaporation of icy grain man- 



tles at temperatures higher than 90-100 K (see e.g. Bottinelli 
|et al.|2004a|b[[Ceccarelli et al.|2007| l, but also species indicative 



of a warm carbon-chain chemistry (CCH and C-C3H2; Sakai 
et al. 2009a b"). Watanabe et al. suggest that 1RS7B could be an 
intermediate-type source between the typical "hot corinos" and 
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"warm carbon-chain chemistry sources". However, a different 
possibility could be that the observed chemical composition is a 
result of photodissociation caused by bright FUV radiation from 
the nearby star R CrA further underscoring the importance of 
understanding the properties of the protostar's environment. 

In this paper we present high angular resolution millimetre 
maps of the line emission in the region surrounding R CrA IRS7 
and aim to quantify the possible physical and chemical impact of 
the iiTadiation by the Herbig Be star R CrA. The paper is laid out 
as follows: Sect. |2] describes the details of the observations and 
reduction techniques. Sect. [3] shows the results of these obser- 
vations. These results are then analysed with LTE and non-LTE 
modelling in Sect. [4] Finally, Sect. [5] discusses the implications 
of the derived physical and chemical properties. 



2. Observations 

2.1. SMA 1 .3 mm observations 



The region around IRS7A, 
SMM lA, and R CrA was 



Submillimeter Amy (SMA; Ho et al. 



IRS7B (SMM IB), SMM IC, 
observed at 1.3 mm using the 
in its compact con- 



2004U 



figuration on 20 August 2006. Six antennae were in the array 
covering baselines from 5-50 k/l. The region was covered by 
two pointings, one located on IRS7B, and one south of R CrA. 
Each of the two pointings covers a circular area with a diame- 
ter of about 7 500 AU dictated by the SMA primary beam size 
of 58", but the two pointings overlap somewhat. The exact co- 
ordinates of the pointings can be found in Table |2] The SMA 
1.3 mm continuum is overplotted on the SCUBA 850 ^im contin- 
uum ( Nutter et al.|2005 ) together with the SMA primary beams 
at the wavelength of the observations in Fig.|2] 

Both pointings were observed with the lower sideband cov- 
ering 216.849-218.831 GHz, and the upper sideband covering 
226.849-228.831 GHz. The SMA correlator was configured to 
give uniform spectral coverage over each of these 2 GHz band- 
widths. The resulting frequency resolution of 0.81 MHz corre- 
sponds to a velocity resolution of 1.1 km s ' at the observed fre- 
quencies. The data were calibrated using the MIR package (| Qi] 
2012[l dedicated to SMA observations. The complex gains were 



calibrated through observations of the quasar 1924-292 with the 
bandpass and absolute levels calibrated through observations of 
the strong quasar 3c454.3 and Uranus. The resulting synthetic 
beam size of the observations was 6.0"x2.9" with a position an- 
gle of 5.3°. 

Further imaging of the data was performed with the MIRIAD 
package ( Sault et aL]|1995 1. Fkst, a linear fit to the continuum 
was subtracted from the (m, v) files. The sky brightness distribu- 
tion was then produced from an inverse Fourier transform of the 
(m, v) data sets, and consecutively CLEANed. Finally, the data 
from the two separate pointings were combined by use of the 
linmos (linear mosaic) task. The 1.3 mm continuum data were 
produced by averaging the (m, v) data across both SMA side- 
bands. The contribution from the strong line emission was esti- 
mated to be less than 10% in the observed continuum peaks (in- 
cluding SMM lA). Previous papers by |Chen & Arce| ( pOTO] l and 
[Peterson et al.| ( [201 \\ have discussed the continuum data from 
these observations separately; here we focus predominantly on 
the line observations. 



' The Submillimeter Array is a joint project between the Smithsonian 
Astrophysical Observatory and the Academia Sinica Institute of 
Astronomy and Astrophysics and is funded by the Smithsonian 
Institution and the Academia Sinica. 
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Fig. 2. SMA 1.3 mm continuum (contours at 3cr, 6cr, 9cr, and 
then increasing in steps of lOcr, where cr - 4.7 mJy beam"') and 
SCUBA 850 |,im continuum (colour scale). The two circles show 
the primary beams of the two SMA pointings. The box shows the 
coverage of the APEX short-spacing data. The synthetic beam 
size of the SMA observations is 6.0"x2.9". Refer to Fig{T]for a 
guide to the symbols used for the compact objects. 

Table 2. Pointing centres and frequency coverage for the SMA 
and APEX observations. 



Telescope 


R.A. 


Dec. 


Line data 


Region 


(12000.0) 


(12000.0) 


[GHz] 




SMA 










IRS7A 


19:01:53.30 


-36:57:21.0 


217.84 


1.0 


IRS7A 


19:01:53.30 


-36:57:21.0 


227.84 


1.0 


IRS7B 


19:01:56.40 


-36:57:27.0 


217.84 


1.0 


IRS7B 


19:01:56.40 


-36:57:27.0 


227.84 


1.0 


APEX 










IRS7B'' 


19:01:56.40 


-36:57:27.0 


217.94 


0.9 



Notes. '"'7x7 pointings in an equatorial grid centred at the tabulated 
coordinate. The separation between the grid points is 14". 



The SMA line data were subsequently combined with APEX 
short-spacing data (see below). 

2.2. APEX 1.3 mm observations 

To study the extended molecular emission in the region, short- 
spacing data for the SMA observations were obtained with 
APEXrl The region was observed in a 7 x 7 pointing pattern 
centred at IRS7B, covering most of the SMA primary beams 
(Fig|2j the westernmost part of the region covered by SMA was 
not covered in the APEX map, since no signal was detected in 
the SMA maps at those positions) in the frequency range be- 
tween 216.951 GHz and 218.849 GHz (corresponding to the 
SMA lower sideband). The beam size of the APEX telescope 



- The Atacama Pathfinder Experiment (APEX) telescope is a col- 
laboration between the Max Planck Institute for Radio Astronomy, the 
European Southern Observatory, and the Onsala Space Observatory 
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Table 1. Continuum point sources in the SMA, Spitzer, and SCUBA data. 



Source name 


RA 


Dec 


Other names" 


Coordinate reference* 




(J2000.0) 


(J2000.0) 






SMM IC 


19:01:55.29 


-36:57:17.0 


SMA 2 (1,2), B9 (3), Brown 9 (4) 


SMA 


IRS7A 


19:01:55.32 


-36:57:21.9 


IRS7W (5), IRS7 (2,6) 


Spitzer 


IRS7B 


19:01:56.40 


-36:57:28.3 


SMM IB (6,7), SMA 1 (2), IRS7E (5) 


SMA 


SMM 2 


19:01:58.54 


-36:57:08.5 


CrA-43 (1), WMB 55 (3,6) 


Spitzer 


RCrA 


19:01:53.67 


-36:57:08.0 




Spitzer 



Notes. References: (1) 



Peterson et al. 



p006) , (6) |Nutter et ar] ( |2605) , (7) |Chen & ArceH2010> . 

{b) Xhp fnorHiniitpc mmp frr\m ^\4A iinrl S^nitv^r rrintiniii 



1 201 l|l (2) [Groppi et al.| j2007| >, (3) |Choi et aL]j2008t , (4) [Miettinen et al.| ( |2008| , (5) [Forbrich etal 



The coordinates come from SMA and Spitzer continuum reported in Peterson et al. 1 201 1 1 



at these frequencies is 28.6", and the individual pointings were 
therefore separated by 14" to give a fully Nyquist-sampled map. 
The frequency resolution of the APEX observations (0.12 MHz) 
corresponds to a velocity resolution of 0.17 km s"^ 

The APEX data were reduced with the GILDAS CLASS 
packag^ Sinusoidal or, in some cases, polynomial (n < 5), 
baselines were carefully calculated and subtracted. The data 
were then exported to FITS format, to allow for import of the 
data to MIRIAD. 

The APEX data were smoothed to match the velocity resolu- 
tion of the SMA data with the regrid task in MIRIAD. The two 
data sets were then combined with the MIRIAD immerge task, 
which linearly merges two sets of data with different spatial reso- 
lution (i.e. applies short-spacing corrections from the single-dish 
data onto the interferometry data). In the combination, primary 
beam correction was applied, which amplifies the noise towards 
the edges of the primary beam. The estimated errors in figures 
in this paper are calculated at the edges of the primary beams. 
However, in the estimates of fluxes in individual positions in the 
maps, the RMS noise in that particular point has been used for 
the error analysis. 

This linear method (immerge in MIRIAD) of combining 
single-dish and interferometry data was chosen in favour of the 
non-linear alternative (mosmem in MIRIAD). The reason for this 
is that the separation between the two SMA pointings is greater 
than 0.5 times the primary beam FWHM, making the mosaic 
somewhat undersampled. The non-linear method would thus risk 
to amplify noise, in particular near the primary beam edges. A 
comparison between the two methods also showed this, with flux 
differences of less than 15% in the inner parts of the primary 
beam, but the non-linear method over-estimating the signal close 
to the primary beam edges. 

2.3. Spitzer, JCMT/SCUBA, and Herschel data 

In addition to the described SMA and APEX observations we 
utilise continuum maps from mid- and far-infrared space-based 
observations from the Spitzer Space Telescope ( [Peterson et al.| 
|20I l| l and Herschel Space Observatory (Lindberg et al., in prep.) 
as well as ground-based submillimetre observations from the 
SCUBA camera at the James Clerk Maxwell Telescope ( |Nutter| 
let al.|2005| l. 



Table 3. Identified molecular Unes in the SMA/ APEX data. 



Species 


Transition 




Frequency" 










[GHz] 


[K] 


SiO 


/ = 5-4 




217.10498 


31.26 


DCN 


/ = 3-2 




217.23854 


20.85 


"CN 


N = 2-1, J = 


3/2-1/2* 


217.26464 


15.66 


"CN 


N = 2-1, J = 


5/2-3/2^ 


217.46715 


15.68 


C-C3H2'' 


6l6-5o5 




(217.82215 


38.61 


C-C3H2'' 


6o6-5l5 




1 217.82215 


38.61 


C-C3H2 


5 14-423 




217.94005 


35.42 


C-C3H2 


524-4l3 




218.16044 


35.42 


H2CO 


3o3-2o2 




218.22219 


20.96 


HC3N 


J = 24-23 




218.32472 


130.98 


CH3OH 


422-3 12 




218.44005 


31.60 


H2CO 


322-221 




218.47563 


68.09 


H2CO 


321-220 




218.76007 


68.11 


'^CN'' 


A' = 2-1, 7 = 


5/2-3/2/ 


226.87478 


11.35 




N = 2-1, / = 


5/2-3/2« 


226.88742 


11.35 




^ = 2-1,7 = 


5/2-3/2'' 


226.89213 


11.35 


C-C3H2'' 


432-321 




227.16913 


29.07 


HC3N" 


J = 25-24 




227.41890 


141.90 



Notes. Frequencies and energies were acquired from the CDMS 
database (Miiller et al. 20011, except for the C-C3H2 data, which are 
from the JPL database (P ickett et al.ri998|^ . 

W Fi = 1-0, F = 0-1 Blend of Fi = 3-2, F = 4-3; F, = 3-2, F = 
3-2; and F, = 3-2, F = 2-1. This is a blend of two C-C3H2 lines. 

Only SMA data are available for these lines. Blend of F = 5/2-3/2, 
F = 7/2-5/2, and F = 3/2-1/2. <«> F = 3/2-3/2. F = 5/2-5/2. 



3. Results 

Eighteen spectral lines from seven molecular species were de- 
tected in the 216.8-218.8 GHz and 226.8-228.8 GHz SMA data. 
The detected species are H2CO, CH3OH, C-C3H2, HC3N, DCN, 
CN (both i^CN and '^CN), and SiO (see Table[3]i. Most of these 
lines, in particular the H2CO and CH3OH lines, show very ex- 
tended emission, which in most cases does not coincide with 
any of the bright compact continuum sources associated with the 
main YSOs, but does align with e.g. the SCUBA peak SMM 1 A, 
south of the suggested YSOs and is also detected in a region be- 
tween R CrA and IRS7|3 

In Fig. [3] the integrated intensity maps of the H2CO 3o3-2o2 
line using SMA data, APEX data, and the combined data set are 



^ GILDAS CLASS (Continuum and Line Analysis Single-dish 
Software) is developed by the IRAM institute, Grenoble, France: 
http : //www . iram . fr/IRAMFR/GILDAS 



Note that Chen & Arce 



1 2010 1 report that no molecular line emis- 



sion from SMM 1 A is seen in the same data; contrasting what is seen in 
the results presented here from our reduction of the data. 
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compared. The bulk of the emission is extended and thus com- 
pletely resolved-out in the pure interferometry data. Combining 
the SMA data with the large-scale APEX maps is therefore criti- 
cal in order to derive quantitative numbers from the data set. All 
12 spectral Hues detected in the SMA 216.8-218.8 GHz band 
were also detected in the APEX data. The SMA and APEX data 
were combined following the methods describe d in|Weifi et al.| 
pool ) and Takakuwa et al.| ( i2003 ) (see also Sect.|2.2[). hi Fig.|4| 
the pure interferometry fluxes and the short-spacing-corrected 
fluxes of the H2CO 3o3-2o2 line from six positions in the inter- 
ferometry map are compared. The pure interferometry lines are 
not only fainter, they also show absorption in several positions, 
since the real fluxes in those positions are lower than the average 
of the extended flux in the field. Typically, less than 25% of the 
peak flux is recovered in the interferometry data. 
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Fig. 4. Comparison between SMA (blue) and SMA/APEX 
(black) H2CO 3o3-2o2 spectra in SMM lA, SMM 2, R CrA, 
IRS7A, SMM IC, and IRS7B. The SMA spectra have been 
scaled by a factor of 4, and there is a 10 Jy beam"' vertical shift 
between the spectra. Units are in Jy beam"' for the SMA beam 
(6.0"x2.9"). 



Figs. |5] and A.1-A.13 show molecular line emission in the 
combined SMA/APEX maps in most of the area covered by 
both the APEX map and the SMA primary beams. Since primary 
beam correction is applied, all contours correspond to the same 
flux level, but not to the same S/N level. Thus, in the channel 



Table 4. Integration intervals and contour levels of the integrated 
intensity maps in Fig. |5] All contour levels are at a 3cr level at 
the edge of the primary beam. 



Spectral line 




V 

[km s"'] 


Contour level 
[Jy beam" 'km s"'] 






r-o ?■ +8 71 


1.46 


H2CO 321-220 




[+3.0; +7.5] 


1.03 


HtCO 322~22i 




[+2.6; +8.2] 


1.15 


CH3OH 422-3 12 




[+2.8; +7.3] 


1.03 


HC3N J = 24-23 




[+4.5; +7.9] 


0.89 


C-C3H2 6l6-5o5, 606 


-5 15 


[+3.7; +7.1] 


0.89 


C-C3H2 5 14-423 




[+4.8; +7.0] 


0.73 


C-C3H2 524-4l3 




[+3.5; +5.7] 


0.73 


^^CNN = 2-l,J = 


3/2-1/2" 


[+4.9; +8.3] 


0.90 


DCN J = 3-2 




[+4.8; +7.0] 


0.73 


SiO 7 = 5-4 




[-2.1; +11.4] 


1.80 


'2CNA/' = 2-l,y = 


5/2-3/2* 


[+2.9; +9.4] 


1.22 



Notes. <°> Fi = 1 

and F = -ill-Ill. 



-0, F = 0-1. <*' Blend of F = 5/2-3/2, F = 7/2-5/2, 



A.l 



A. 13 I the 0.5 Jy beam contours are approxi- 



maps (Figs.^ 

mately at a 3cr level on the primary beam edges (see Fig.|2]l, but 
at a 6cr level in the beam centres. 

The molecular lines seen in the SMA data can roughly be 
separated into three groups by their distribution in the R CrA re- 
gion. The first group consists of all detected H2CO and CH3OH 
lines. The second group consists of the HC3N, C-C3H2, DCN, 
'^CN, and '^CN lines. Finally, the distribution of the SiO emis- 
sion is not similar to those of the other molecular species. 

3.1. HiCOandCHiOH 

In the APEX data, the H2CO 3o3-2o2 line has a non-Gaussian 
shape not seen in the higher-excited H2CO lines, suggesting ei- 
ther that the emission has multiple components, or that optical 
depth effects are present (see Fig.[6]l. 

The H2CO and CH3OH lines show strong emission in two 
ridges; one south of IRS7B, with one of the peaks in this ridge 
almost coinciding with the suggested starless core SMM lA 
dChen & Arce 2010); the other between SMM IC and R CrA 
(see Figs. |5] and A.l -A.4). Both these regions are relatively ex- 
tended also in the velocity domain (vlsr = 4 to 7 km s"' for the 
northern ridge and vlsr = 5 to 8 km s"' for the southern), but 
with individual clumps that peak at different velocities in these 
intervals. Note, however, that significant H2CO 3o3-2o2 emis- 
sion is seen all over the field between 4 and 7 km s"' . This mor- 
phology is different from the sources studied by |j0rgensen et al. 
( |2007[ ), where H2CO and CH3OH emission mostly is found to 
be associated with the YSOs and their outflows. 

The H2CO 3o3-2o2 emission is of comparable strength in the 
northern and southern ridges, whereas the higher excited H2CO 
lines and the CH3OH line are stronger in the southern ridge. The 
difference in the H2CO line intensities suggests a higher temper- 
ature in the southern ridge than in the northern or alternatively 
a higher optical depth of the main transition. The emission peak 
in the southern ridge does not coincide with the centre of the 
extended SCUBA dust emission, which lies approximately 10" 
southwest of the strongest H2CO 3o3-2o2 peak in the southern 
ridge. However, the SCUBA dust emission peak coincides well 
with a secondary peak of the higher excited lines (see Fig.|7]). 
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Fig. 3. Left: H2CO 3o3-2o2 line emission detected by the SMA integrated between and +9 km s"'. Middle: H2CO 3o3-2o2 line 
emission observed by the APEX telescope in a 7x7 grid integrated between and +9 km s"' . Right: H2CO 3o3-2o2 line emission in 
the combined SMA interferometry and APEX short-spacing data integrated between and +9 km s 

Contours are at 2.5 Jy beam"' km s"' levels (converted to SMA beams, ~ 5cr at the edge of the SMA primary beam) in all three 
maps. Negative intensities are represented by dashed contours. Coordinates are offsets from the central position of IRS7B. Refer to 
Fig|T]for a guide to the symbols used for the compact objects. 



For the lower LSR velocities (-1 to +3 km s most H2CO, 
CH3OH, and SiO emission is seen in the northeastern part of 
the field covered by the SMA/APEX observations. This low- 
velocity emission is not connected to the emission around IRS7 
- if that was the case, a much smoother westward transla- 
tion of the emission with increasing velocity would be seen in 
Fig. A.l Instead, the H2CO emission shows a distinct jump 



around 3 km s . This indicates that the northeastern emission 
originates in another source, outside or on the edge of the field- 
of-view, perhaps SMM 2. This emission is also fainter than that 
in the centre of the field. 



3.2. HC3N, C-C3H2, DCN, and CN 

The second group of lines (HC3N, C-C3H2, DCN, and CN) 
shows relatively weak emission. Except for CN, they are all 
present in an EW-extended region between IRS7B and SMM 1 A 
(see Figs. [5] and 



A.5 - A.12i. C-C3H2 and CN are the only 



molecules in the SMA/APEX data that have their strongest emis- 
sion associated with the YSOs: C-C3H2 peaks at IRS7B, and ex- 
tends in a northward collimated outflow, whereas '^CN shows 
a strong clump just north of IRS7B. Extended HC3N emission 
is detected south of IRS7A and southwest of IRS7B. The '^CN 
emission is rather diffuse and extended, but also of quite low 
S/N. For the '^CN lines, no APEX short-spacings are available, 
so no conclusions about the extended emission can be drawn 
from this map. However, a strong clump of '^CN emission is 
seen peaking just north of IRS7B, and a fainter clump is seen 
between IRS7A and R CrA. 



SMM 2 at a velocity slightly lower than the systemic velocity of 
R CrA (0 to 3 km s"'). The latter two emission regions are out- 
side the SMA primary beams, and thus only tentative. All SiO 
peaks coincide with relatively faint H2CO emission, which is in 
agreement with SMA observations of other prominent embed- 
ded protostars showing these molecules to coincide in outflow- 
affected gas ( J 0rgensen et al.|2007| l. 

The most interesting SiO feature is the point source south of 
R CrA, since it together with R CrA falls on a straight line with 
several Herbig-Haro (HH) objects. Hartigan & Graham (1987|l 
observed several HH objects (HH98, HH98B, HHIOO, HH97, 
HH96, and HHlOl) in a straight line that extends southwest- 
wards from R CrA, through the SiO peak in the SMA data. 
They argue that all HH objects are driven by the infrared source 
IRSl, almost 1' southwest of R CrA. However, Wang et al. 
(2004) argue that, considering the positions of the objects, the 
outflow might as well be driven by R CrA. The proper-motion 
velocities of the objects southwest of IRSl were measured by 
[Peterson et al.| ( [20II| l, and are indeed directed southwestwards. 
The proper-motion velocity of HH98B (MHO 200 IQ) is also 
measured, which is directed westwards, and suggested to be 
driven by SMM IC (SMA 2), and is thus not related to the other 
HH objects along the line. 

As a summary, the Herbig-Haro objects in the region have 
several possible driving sources (including IRSl, SMM IC, and 
R CrA), and the SiO peak south of R CrA could be associated 
with any of these objects. There are no sources described in 
that can explain the SiO emission asso- 



Peterson et al. 



(2011 



ciated with SMM I A, nor have any HH outflows been observed 
in its vicinity. 



3.3. SiO 

The SiO emission is present in four distinct regions separated 
both in position and velocity (see Fig. A.13i. It appears as a 
bright point source south of R CrA, west of SMM IC, at ve- 
locities (9 to 12 km s ') considerably higher than that of R CrA 
(5 km s"'); in somewhat extended emission around SMM lA at 
vlsr = 6 km s ' ; as a pointlike source west of T CrA at low ve- 
locity (-1 to +2 km s"'); and finally some faint emission east of 



4. Analysis 

With the high resolution interferometiic maps from the SMA and 
the recovery of the more extended structure through the APEX 
observations, it is possible to investigate the physical conditions 
in the region around the protostars close to R CrA. In the follow- 
ing we constrain the physical properties in the region and subse- 
quently (Sect. |5]l the implications concerning the radiation field 



6 



J. E. Lindberg and J. K. J0rgensen: Strong irradiation of protostellar cores in CrA 




30 20 10 -10 -20 -30 -40 30 20 10 -10 -20 -30 -40 30 20 10 -10 -20 -30 -40 



R,A offset (") 

Fig. 5. SMA/APEX integrated intensity line maps of all spectral lines in the SMA coverage except for '■'CN = 2-1, J - Sjl-'ill 
(blend of several hyperfine transitions), which is barely detected in the interferometry data. The plotted '^^CN line has the quantum 
numbers Fx - 1-0, F - 0-1. No APEX short-spacing data were used for the '^CN - 2-1, J - 5/2-3/2 (blend of several hyperfine 
transitions) line. All contour levels are at a 3cr level at the edge of the primary beam. Refer to Table|4]for the exact contour levels 
and the integration intervals. Refer to Fig{T]for a guide to the symbols used for the compact objects. 

around the embedded protostars and in particular the influence 
from R CrA itself. 
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Fig. 7. Integrated intensity of the H2CO 3o3-2()2 line (left) and H2CO 322-221 line (right). Contours every 2 K km s"' (~ 3cr). Both 
lines are integrated between and 9 km s ' . Greyscale background of both figures is SCUBA 850 ^m continuum. Refer to Fig{l]for 
a guide to the symbols used for the compact objects. 



4.1. LTE estimates of physical properties 

Observations of several rotational transition lines of the same 
molecule can be used to construct a so-called rotational dia- 
gram for a given molecular species from which the column den- 
sity and rotational temperature of the species can be estimated 
( [Goldsmith & Langer|1999) . If one assumes that the molecules 
are excited under local thermal equilibrium (LTE) conditions and 
that the lines are optically thin, the rotational temperature is an- 
other measure of the kinetic temperature of the species. In the 
case of H2CO, the three lines at 218 GHz make it possible to cal- 
culate first-order approximations of the temperature and column 
density (see e.g. |Mangum & Wootten|1993[|Maret et al.|2004| l. 
However, if the lines are optically thick, the rotational diagram 
method will usually overestimate the kinetic temperatures, since 
the lower-excited lines will be more severely affected by the op- 
tical thickness than the higher-excited lines. 

The temperature and column density of the H2CO in the re- 
gion covered by the SMA and APEX data is estimated by the 
use of these three spectral lines in all parts of the SMA primary 
beams where the integrated intensities between and 9 km s"' 
of all three spectral lines have a signal-to-noise ratio of at least 5. 
In Fig. [8] the rotational temperatures and column densities esti- 
mated from the SMA/APEX maps of the H2CO 3o3-2o2, 32i-22o, 
and 322-221 lines are shown. Note that the lines in this study are 
para-H2CO lines, and to find the total H2CO column densities, 
the ortho/para ratio must be accounted for. In this paper, the p- 
H2CO column density will be stated. Assuming an ortho/para 
ratio of 1.6 (Dick ens & Irvine|19 99; J0rgensen et al. 2005]l, the 
total H2CO column density can be estimated by multiplying the 
P-H2CO column density by 2.6. 

The maps show a wide spread in rotational temperatures 
throughout the region; from 30 K in the colder regions (mostly 
associated with low H2CO intensities, and thus low S/N lev- 
els), up to 120 K in a few peaks close to SMM lA. From the 
SMA/APEX map of the IRS7 cloud (Fig.|8]l, it can be found that 
the highest rotational temperatures are not associated with the 
YSOs, but with the southern ridge around SMM lA, and also 
with the northern ridge close to R CrA. In some points in the 
southern ridge, a rotational temperature of 120 + 20 K is found. 



Using the APEX pointing centred at IRS7B, a single-dish ro- 
tational temperature of 62 K was calculated in a 28.6" beam 
around IRS7B. However, the highest temperatures are likely op- 
tical depth effects, and the physical temperatures should rarely 
exceed 70 K (see Sect. |4.1.1|). The estimated column density 



varies between 10'^ and 10'"* cm"^ in the region, and is fairly 
well-correlated with the rotational temperature. As can be seen 
in Fig. |8] the rotational temperature is not peaking at the same 
place as the H2CO 3()3-2o2 (main line) emission. The rotation 
temperatures at the positions of the YSOs are relatively moder- 
ate, around 40-50 K (see Fig.|9]l. 

Only one of the three high-temperature H2CO peaks around 
SMM lA coincides with one of the three condensations within 
SMM lA suggested by [Chen & Arcel ( [20T0) l: SMM lA-a, the 
weakest 1.3 mm continuum source within SMM lA, coincides 
with the peak with the second highest rotational temperature. 
The much more luminous SMM lA-c is found in a point with 
relatively low H2CO temperature (T x 65 K). The northern 
ridge of extended H2CO emission has a more even and some- 
what lower rotational temperature {T ^ 40-50 K) than the south- 
ern ridge (r 50-120 K). 

These temperature and column density estimates suffer from 
four main sources of eiTors: The statistical eiTor (noise); the cal- 
ibration error of the telescopes; the primary beam error, which is 
notable far from the primary beam centres; and the errors arising 
from the model assumption of LTE conditions. The statistical 
uncertainty, including the primary beam error, is estimated from 
the RMS of the noise in the combined SMA and APEX data. The 
systematic calibration error is estimated to 20% of the line flux 
(note that both SMA and APEX will introduce systematic errors, 
which both should be in the order of 20%, but since the line flux 
is dominated by the APEX data (see Fig. |4]), it can be assumed 
that the total calibration error is almost equal to the contribu- 
tion from APEX). However, this error does only apply to the 
column density estimates, since the rotational temperatures are 
proportional to the logarithm of the line ratios and thus only de- 
pend on the relative errors in the line flux estimates due to noise. 
Typically, the rotational temperatures are therefore determined 
to about 10% accuracy given the model assumptions. 
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Fig. 8. Left: H2CO rotational temperature map of the IRS7 region from SMA/APEX data. Temperatures are in K. Right: P-H2CO 
LTE column density map of the IRS7 region. The colour scale bar shows values in logjg A^, where A^ is given in cm"^. 
Both maps use the integrated intensities of three para-H2CO lines around 218 GHz, and only data points that are inside the SMA 
primary beams and where all three lines have 5cr detections or stronger are considered. All other points are represented by white 
pixels. The temperatures and column densities are calculated with rotational diagrams for line intensities integrated between and 
9 km s"'. The contours show the integrated intensity of the H2CO 3()3-2o2 line in the same velocity interval, with contours every 
2 K km s ' (~ 3cr at the edge of the primary beam). Refer to Figlllfor a guide to the symbols used for the compact objects. 



|Maret etaT] (|2004') studied the H2CO abundances in a sam- 
ple of Class protostars using the IRAM 30 m telescope and 
the JCMT. Several P-H2CO and 0-H2CO spectral lines between 
140 GHz and 365 GHz were observed in eight sources, and the 
highest H2CO rotational temperatures were found in three hot 
corino sources in NGC 1333 (25^0 K), whereas the rotational 
temperatures in the other sources were in the order of 10-20 K. 
However, if only using the 218 GHz lines, much higher rota- 
tional temperatures (up to several 100 K) are reached in the 
hot corino sources, but with large uncertainties. Furthermore, 



J0rgensen et al.| ( j2005| 2007) demonstrated that the single-dish 
H2CO observations are difficult to interpret for the specific 
sources - in particular, that the high rotational temperatures are 
a likely result of extended outflow emission picked up by the 
single-dish observations. 



4.1.1. Optical depth 



To investigate whether the assumption of optically thin H2CO 
lines is valid, the intensities of the Hj^CO 3i2-2n line at 
225.6978 GHz and the H'^CO 3i2-2ii line at 219.9085 GHz 
(both from a large APEX survey; Lindberg et al. in prep.) were 
compared. The ratio between the integrated intensities of the two 
lines was calculated to be 35 (Fig 

pared to the typica l ISM '^C/'^C ratioTn the local ISM of 77 ± 7 
( [Wilson & Rood|[l994) . The ratio measured in IRS7B is thus 
a factor of 2.2 lower, giving an optical depth of the o-Hj^CO 
3i2-2ii line of t ^ 0.8. Using equations for the total column 
density of a molecule (see e.g. Goldsmit h & Langer|1999), and 
assuming an ortho/para ratio of 1.6 ( |Dickens & Irvine 1999[ 
J0rgensen et al.|2()05 1, and Tex = 62 K (the rotational tempera- 
ture observed in the APEX beam), we find that this corresponds 
to T * 0.69 for the P-H2CO 3o3-2o2, and r a; 0. 18 for the P-H2CO 
321-220 and P-H2CO 322-221 lines. Note, however, that these op- 
tical depths are averages over the APEX 28.6" beam, and that 
locally elevated optical depths are likely. 

Since the optical depths of the three spectral lines used for 
the rotational temperature calculation are different, the rotational 
temperature will be overestimated. When compensating for op- 



tical depth effects using the values estimated above and the opti- 
cal depth correction factor given by Goldsmith & Langer ( 1999) , 
the rotational temperature in the APEX beam drops from 62 K 
to 47 K. Since the distribution of the optical depth across the 
SMA/APEX map is unknown, non-LTE methods must be used 
to acquire a more accurate temperature estimate, but applying 
the large-scale optical depth on the calculations of the rotational 
temperature peaks in the southern ridge makes these drop from 
100-120 K to 65-75 K. 



4.2. Non-LTE analysis 

For further analysis of the emission, we utilised the RADEX es- 
cape probability code ( van der Tak et al.|2007 1 to solve the exci- 
tation of H2CO under non-LTE conditions assuming a uniform 
spherical molecular cloud. Besides allowing for a more detailed 
comparison of the importance of the density and temperature, 
this method also provides a path to take optical depth effects 
into account. Collision rates come from Green ( 1991 1, and were 
retrieved from the LAMDA database ( |Schoier et al.|2005| l. 

In the absence of a central heating source, the temperature 
of a molecular cloud can hardly exceed 10-15 K assuming no 



external heat contributi on, and typically the density would be 



10 1, which should be com- ~ 10 -10 cm (Bergin & TafaUa|2007 |. Furthermore, the low 

mass YSOs should not be able to heat the gas on scales observ- 
able in the SMA/APEX map. Consequently, assuming external 
heating, a fairly even temperature distribution is expected across 
each dense ridge. Thus, it is likely that the wide spread of tem- 
peratures found in the rotational temperature map (Fig. |8]l is 
caused exclusively by column density effects, and that each of 
the two ridges have uniform temperatures. To exclude all emis- 
sion outside these ridges, only pixels where all three H2CO lines 
are detected to at least an 8cr level are considered in the RADEX 
calculations. 

In the ridges, the H2CO 321-220/H2CO 322-221 ratio has an 
average of 0.98, and the ratio lies between 0.8 and 1.2 in more 
than 90% of the pixels in the two ridges. The average optical 
depth of the H2CO 3o3-2o2 line in an APEX beam around IRS7B 
is 0.7 (see Sect. |4.1.1| i, but might be several times higher in the 
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Fig. 6. APEX spectra for H2CO 3o3-2o2 (black) and H2CO 
322-221 (blue), the latter line is scaled with different factors (4.0, 
4.7, and 5.1 from top to bottom) to make the Gaussian tails 
fit. The vertical shifts between the spectra are 5 K. The H2CO 
321-220 (not displayed) has a shape similar to the H2CO 322-221 
line. 



peaks. These parameters are found to be relatively temperature- 
insensitive - at 30-90 K they require « ^ 10^ cm"^ and ~ 
lO'^-lO'^ cm-2. Where the H2CO 321-220/H2CO 322-221 ratio is 
greater than 1, « <: 10^ cm"^. This suggests that the H2 density in 
both ridges must be relatively high compared to what one would 
expect in a prestellar cloud. 

RADEX calculations are run on each of the individual pixels 
in the ridges to investigate which uniform temperature fits the 
data best. This non-LTE model is tested with uniform H2 den- 
sities of 10^ 10^ 10\ and 10** cm , and the column density 
is set to be a free parameter across the cloud. The two ridges 
are treated separately, and the most likely temperature of each 
ridge is established by calculating the of the intensities mod- 
elled by RADEX compared to the observed intensities, and then 
minimising the median of these values for all pixels in each 
ridge. The temperatures where the optical depth of any of the 
lines is negative are excluded (a negative t corresponds to mas- 
ing, which is very unlikely to occur on these large scales, which 
removes some high-temperature solutions). Temperatures where 
the optical depth of the H2CO 3()3-2o2 line is greater than 10 are 
also discarded (this is more than an order of magnitude larger 
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Fig. 9. H2CO rotation diagrams for the YSOs in the SMA/APEX 
field (SMM IC, IRS7B, 1RS7A and SMM lA). 
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Fig. 10. H'^CO 3i2-2ii (black) and H'^CO 3i2-2ii (blue, scaled 
by the factor 35) in 1RS7B. Data from APEX line survey of 
IRS7B (Lindberg et al., in prep.). 



4.1.1 which excludes 



than the average value estimated in Sect, 
some low-temperature solutions). 

From the RADEX calculations, it is found that the density 
must be higher than 10^ cm"^, since densities lower than that 
require a very high optical thickness (r > 10) to produce the 
observed intensities in several positions at reasonable temper- 
atures (T < 90 K). The optical depths found at 10* cm"^ are 
more likely, but for the highest temperatures (T ^ 70 K), one 
or several of the lines would be masing. The lowest-^^ temper- 
atures (excluding r < and t > 10) for each of the fits can 
be found in Table |5] where it can be seen that the best fits give 
temperatures of at least 50 K in the southern ridge in any sce- 
nario (see also Fig. [TT) . From Fig. [12] it can also be found that 
given the most likely column densities estimated from the opti- 
cal depth (see Sect. 



4.1.1 



the temperature is higher than 25 K 
in the northern ridge and higher than 30 K in the southern ridge 
with 3cr certainty. Consequently, both the temperatures and den- 
sities in the region are likely considerably higher tha n the 10- 



15 K and 10"^-10^ cm ^ expected in pre-stellar cores (Bergin & 
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|Tafalla|2007| |Ward-Thompson et al.|2007| l. The best-fit colutnti 
densities are shown in Fig. 13 and lie around lO'^-lO'"^ cm"^ 
in both ridges (shghtly higher in the southern ridge than in the 
northern). 



Table 5. Temperatures calculated by the non-LTE code RADEX. 
The given temperature coiTesponds to the lowest-;^'^ fit, and the 
limits correspond to values where < 2.3. The column densi- 



ties are shown in Fig. 13 



n 


r in N ridge 


r in S ridge 


[cm-3] 
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10^ 
10' 
10« 


65%" 

4.7+45 
^'-16 


76+' * 
67+" * 

"'-19 

61+'' 

-2^ 



Notes. r > 10 for the HiCO 3o3-2o2 line for all T where median 
< 2.3. 

**' The higher temperature limit corresponds to the limit where at least 
one of the lines starts to mase (t < 0). 



In all the fits, a dip in the fit around 15 K is found. This 
is not the minimum-;^^^, but still within the Icr limit. However, 
those low-temperature fits always require very high optical 
depths (r > 100 for n < 10^ cm-\ t > 25 for n > 10^ cm^^), 
which is one to two orders of magnitude higher than the opti- 
cal depth averaged over an APEX beam centred at 1RS7B (see 
Sect. |4.1.1|l. If we assume a n H2CO/H2 ratio in the order of 
10"^ ( J0rgensen et al.p 2005 ), and that the source is spherical 
and has a size in the order of 1", the observed line strengths 
would require an H2 density in the spherical cloud in the order 



of 10^-10'" cm assuming the low-temperature solution. The 
high-temperature solution is thus much more likely. 
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Fig. ll.x^ fits to RADEX data and the observed H2CO lines to a 
constant temperature with variable column density (see Fig.[T3]l 
and an H2 density of 10*' cm"^ for the northern ridge (dashed) 
and the southern lidge (solid). The greyed areas in the plot cor- 
respond to solutions with r > 10 (left) and r < (right). 
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Fig. 13. P-H2CO column densities for the least-;^^^ fits to RADEX 
data and the observed H2CO lines to a constant temperature in 
each ridge (60 K in northern lidge and 76 K in southern ridge) 
and an H2 density of lO*" cm"^. The colour bar shows values in 
log JO A^, where is given in cm" . Refer to Fig [I] for a guide to 
the symbols used for the compact objects. 



5. Discussion 

5.1. Internal vs. external heating 

It was not possible to fit the observed line intensities with tem- 
peratures lower than 35 K in either the LTE models or the non- 
LTE models of the physical properties of the cloud. More likely, 
the temperatures in both ridges lie around 50-60 K, and per- 
haps as high as 90-100 K in a few clumps south of SMM lA. 
SMM IC lies on the southern border of the northern ridge, but 
IRS7A and IRS7B do not belong to any of the ridges. In the 
regions between the ridges, the column densities are approxi- 
mately an order of magnitude lower than in the ridges. The tem- 
peratures around the YSOs are likely also somewhat lower (30- 
50 K) than in the ridges, but these estimates are not as reliable, 
since the S/N levels are lower in these regions. 

According to dust radiative transfer calculations, the temper- 
ature of the envelope heated only by an internal 20 L0 Class O/I 
YSO will at most reach at most 30 K on the scales of 1 000 AU 



(see, e.g., J0rgensen et al. 2006 and references therein). For pro- 
tostars with even lower luminosities, such as the YSOs in IRS7, 
the internal heating should be of even lesser importance, and on 
the scales of the resolution of the SMA data (400-800 AU) tem- 
peratures in the order of 50-60 K as a result of internal heating 
are unrealistic. However, if heating from an external UV field is 
applied, the models discussed below show that it is possible to 
reach the measured temperatures. 

Not only does the region south of IRS7 show high H2CO 
temperatures, the 850 (,im dust emission also has a strong peak 
southwest of the strongest H2CO peak (see Fig. [7]). This dust 
emission peak coincides well with the temperature peak (T « 
100 K) in the rotational temperature map (Fig. |8]l. Thus, the 
higher rotational temperature measured in the peaks could also 
be an optical depth effect caused by the high column density. 

The angular separation between R CrA and SMM lA is 
approximately 50", which corresponds to a projected separa- 
tion of 6500 AU assuming the distance of 130 pc to the cloud 
( Neuhauser & Forbrich|2008). If the lumin osity of R CrA is as- 
sumed to be ~ 200 L0 ( Bibo et al.|1992| l, its surface tempera- 
ture should be in the order of 10 000 K. Assuming a black body 
with this temperature, slightly less than 50% of its luminosity 
(~ 90 Lq) would fall in the UV range between 300 nm and 
13.6 eV. Such a star would increase the interstellar radiation field 
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Northern ridge, = 5x10^'' cm 



Southern ridge, A'' = 5x10^^ cm 




Fig. 12. -confidence contour plots for temperature and number density determined from the observed H2CO Unes and RADEX 
models. Left: Fit to the average intensity in the northem ridge, N = 5x lO'-' cm"^. Right: Fit to the average intensity in the southem 
ridge, N = 5 x 10'^ cm"^. The greyed areas represent solutions with masing H2CO emission. 



(ISRF) by a factor of ;ifisRF ~ 3 000 at the distance 6 500 AU, as- 
suming an interstellar radiation field of 8 x 10""* erg cm"^ s"' 
(twice the Habing flux). For a slightly lower luminosity for 
R CrA (100 Lq is the lower hmit given by i Bibo et al.J1992[ ), and 
the distance being somewhat higher (if R CrA and SMM 1 A are 
not in the same plane of sky), xisrf ~ 750. Such high UV fluxes 
would create a Photon Dominated Region (PDR), which would 
give rise to high CN abundances, in agreement with Watanabe 
|etaL| ( [20T2l i. 

To estimate the relative importance between internal and ex- 
ternal heating, we constructed a 1 -dimensional model for the 
envelope around IRS7B using the dust radiative transfer code 
Transphere ( jPuUemond et al. 2002| l. We assumed a power-law 
density profile for the envelope with n oc /? '-^ and then varied 
the density profile normalisation, envelope size, central source 
luminosity, and interstellar radiation field. We adopted the dust 
opacities from [Ossenkopf & Henning (1994i corresponding to 
coagulated dust grains with thin ice mantles (OH5). Transphere 
calculates the self-consistent temperature in the spherical enve- 
lope with these parameters. For comparison to the observations 
we subsequently used the SKY ray-tracer from the RATRAN 
radiative transfer code (Hogerheijde & van der Tak 2000) to cal- 
culate continuum maps at different wavelengths. To allow for 
comparison with SCUBA and Herschel SED fluxes, the enve- 
lope dust emission at the corresponding telescope primary beam 
edges was subtracted from the peak emission. This difference 
was then compared to measured SCUBA (Nutter et al. 2005| l 
and Herschel (Lindberg et al., in prep.) continuum fluxes. 

The SED produced by the best fit (with both internal and 
external heating) is compared to the SCUBA and Herschel con- 
tinuum fluxes in Fig. [14] The eight wavelengths at which the 
Herschel fluxes were estimated were chosen to cover each of the 
four PACS bands with two data points equidistant from the band 
edges and each other. The systematic continuum data errors are 
estimated to 20% of the fluxes for the Herschel data (Green et 
al., in prep.), and 50% and 20% for the 450 [xm and 850 ^im 
SCUBA data, respectively ( jDi Francesco et al.|2008i l. The best 
fit was produced for a central source of 4 Lq, with an envelope 
between R,^ < 200 AU and R^m = 10000 AU with n ~ ' ^ 
no = 1-3 X 10^ cm"^ at the reference radius Rq - 1 000 AU (cor- 



responding to an envelope mass of 2.2 Mq), and an interstellar 
radiation field with;i'isRF = 750. 



The calculated temperature profile is shown in Fig. 15 where 
r 25 K for * 2 500 AU, the distance between 1RS7B and 
SMM lA. Also shown is the model without the external radia- 
tion field: as argued above, in this model the temperature drops 
below 20 K at radii larger than about 700 AU, thus not reproduc- 
ing the large-scale temperatures observed in SMA/APEX data. 
In contrast, in the model with the enhanced external radiation, 
the temperature never drops below 25 K, and increases reaching 
50 K in the outermost regions of the envelope. With an enhanced 
external radiation field O^isrf = 750) and no internal source of 
heating, the temperature stays above 20 K in the whole envelope. 

Naturally, this 1 -dimensional model is a simplified descrip- 
tion of the physical properties in the region. For example, it is 
unlikely that the density profile continues to follow the power- 
law at large radii. Since we see H2CO emission in the whole 
field, the density should at least exceed the effective density of 
the H7CO tra nsitions, i.e. at least be in the order of 10^ cm"^ 



(Evans 1999 1. Adding an ambient density of 10 cm to the 
density profile only marginally changes the temperature profile 
(Fig.[B). 

The 1 -dimensional model is also a simplification since 
R CrA in reality would cause a radiation field that decreases from 
northwest to southeast, instead of acting constantly from all di- 
rections. The fact that the northern ridge has a lower measured 
temperature than the southern ridge, despite having a smaller an- 
gular separation with R CrA, is most likely a geometric effect. 
Either it is caused by extinction and local density distribution 
differences, or the ridges simply are positioned at different dis- 
tances from us, making the actual distance between the northem 
ridge and R CrA larger than that between the southern ridge and 
R CrA. The difference in LSR velocities between the ridges do 
suggest that they are situated in different planes of sky. Still, this 
simplified model is useful for a first interpretation of the data at 
hand. It shows that the relatively high temperatures in the ridges 
cannot be explained by the intemal heating from the protostars, 
but that they instead can be explained by R CrA acting as an 
external heating source upon the envelope. 
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Fig. 14. The SED calculated from the envelope with a 4 Lq cen- 
tral source and an interstellar radiation field with ^isrf — 

750 

by the use of RATRAN (the solid line in Fig. 15 i compared with 
the real data points from Herschel and SCUBA continuum mea- 
surements. 




Fig. 15. Temperature as a function of radius for the four differ- 
ent models produced with Transphere. Black solid line: With a 
central source of 4 Lq and an interstellar radiation field with 
;fiSRF = 750. Dotted line: With a central source of 4 Lq and 
a normal interstellar radiation field (^isrf - !)■ Dashed line: 
Without a central source, but with an interstellar radiation field 
with xiSY^F - 750. Grey solid line: As the black solid line, but 
with an ambient density of 10^ cm"^. 



5.2. Implications 



As demonstrated 

(Sect, 



5.1 



based on the radiative transfer models 
), the temperature in the protostellar envelope of 1RS7B 
IS cuiTently above 20 K due to the irradiation by R CrA. This may 
have important implications for the chemistry of the envelope 
and thus potentially also the composition of the the emerging 
circumstellar disc. 

The gas-phase chemistry of the protostellar envelope would 
be dominated by radicals such as CN and C2H prominent in 
PDRs or regions of protostellar envelopes affected by strong UV 
irradiation (e.g., Jansen et al. 1995; S ternberg & Dalgarno|1995 



Fuente et al. 



1995, ,j0rgensenn2OO4j . In fact, in the line survey 



Watanabe et al. (2012 1, an over-abundance of these 



oflRSTBby 

species is seen. 

The strong UV irradiation may, however, also affect the 
grain surface chemistry, potentially leading to the formation of 
more complex organic species. If R CrA has been as lumi- 



nous throughout the evolution of IRS7B (from the time of its 
dense prestellar core stage) this would have hindered significant 
freeze-out of molecules such as CO, CH4, and N2, being some 
of the starting points for grain-surface reactions leading to ze- 
roth generation species such as CH3OH ( [Herbst & van Dishoeck 
2009 ). However, if R CrA itself has undergone significant evo- 
lution in terms of luminosity it may just as well have stimulated 
the formation of complex organic molecules throug h the pho- 
tochemistry in CHaOH-rich ices ( |6berg et a"L||2009 1. With the 
present observational data, it is difficult to distinguish these sce- 
narios: the low internal luminosity of 1RS7B means that the re- 
gions where complex organic species will return to the gas-phase 
are small (:$ 100 AU) and severely diluted in typical single-dish 
beams. Future high angular resolution observations with ALMA 
will, however, reveal emission of complex organic molecules if 
present. 



6. Conclusions 

We have presented high angular resolution observations of 
young stars in the region of R CrA from the Submillimeter 
AiTay augmented by short-spacing maps from the APEX single- 
dish telescope. Together the observations reveal the physical and 
chemical structure of the environment of the deeply embedded 
protostars on scales from 400 to 8 000 AU. The main conclusions 
of the paper are: 

1. In the combined SMA/APEX data, we find two extended 
ridges of luminous H2CO emission. They are not associated 
with the YSOs previously observed in the 1RS7 region. The 
H2CO emission extends across the whole field-of-view of the 
observations, spanning more than 10000 AU. 

2. The physical temperatures in the two ridges around 1RS7 
were estimated with both LTE and non-LTE methods. They 
were in both cases found to be no lower than 30 K, more 
likely 50-60 K. These temperatures are too high to be at- 
tributed to internal heating from the low-mass YSOs in the 
region. The peculiar shape of the ridges, and their location 
with respect to the YSOs, also suggest that this region is 
strongly affected by external heating. 

3. The primary suspect for the external heating of the ridges is 
R CrA, a 3 Mq Herbig Be star in the northwestern parts of 
the cloud. The morphology of the H2CO-luminous ridges is 
relatively symmetric around the YSOs, which could indicate 
that the ridges are parts of a spherical shell around the low- 
mass YSOs. 

4. It was possible through radiative transfer modelling to ex- 
plain the elevated temperatures by the UV radiative contri- 
bution from R CrA. As discussed in Sect. [T] earlier studies 
have shown that the evolution of a protostar can be affected 
by high-mass stars in their environments, but this study sug- 
gests that also intermediate-mass stars can have a noticeable 
effect on the physical properties in a whole protostellar enve- 
lope. The model also suggests an internal luminosity of 4 Lq 
and an envelope mass of 2.2 Mq for R CrA IRS7B. 

5. The higher temperature of the southern ridge as compared to 
that of the northern ridge suggests that the southern ridge is 
actually closer to R CrA than the northern ridge is, as a pro- 
jection effect of the ridges being situated in different planes 
of sky. 

6. The strong UV flux from R CrA obviously has an im- 
pact on the gas-phase chemistry on large scales in the re- 
gion. Potentially it could also be critical for the ice chem- 
istry, either inhibiting CO freeze-out or promoting formation 
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of complex organics. The effect on the chemistry depends 
on how early in the evolution of the protostellar envelope 
around IRS7B the luminosity from R CrA started to domi- 
nate the region. The elevated temperature in the cloud could 
also explain the unusual chemistry in the region, since much 
of the CO would evaporate from the dust grains at a very 
early stage, which would make the formation of complex or- 
ganics on the dust grains inhibited. 

This paper demonstrates the importance of understanding the 
nearby environment in the studies of protostars, and also the 
need for high-resolution observations. Unresolved single-dish 
observations of regions similar to IRS7 would not have shown 
the large-structure high-temperature gas found in this study, and 
without the knowledge of these structures, a high H2CO tem- 
perature would likely have been attributed to the central source. 
Single-dish observations might still be useful in the comparison 
of the excitation of different molecules given the prior knowl- 
edge from this work. Thus, we suggest observations of more 
molecular lines to establish whether the high rotational temper- 
ature of H2CO is reflected also in other species. 

With the advent of ALMA, even more high-resolved studies 
of protostars will be made possible, making investigations where 
heating of the envelope from the central source and external stars 
are possible to distinguish. To study the more central and warmer 
parts of the protostellar envelopes, observations of higher fre- 
quency are needed. Such observations are currently being pro- 
vided by Herschel Space Observatory, and will be discussed in 
a future paper (Lindberg et al., in prep.). These observations will 
make it possible to study whether the high level of excitation of 
the molecular gas is mainly caused by a PDR or by shocks. 
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Appendix A: Spectral line channel maps 
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Fig.A.l. SMA/APEX channel maps of the H2CO 3o3-2o2 hne. Contours every 0.5 Jy beam ' km s ' (~ 3cr at edge of primary 
beam). Refer to Fig|T]for a guide to the symbols used for the compact objects. 
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Fig.A.2. SMA/APEX channel maps of the H2CO 32i-22o line. Contours every 0.5 Jy beam ' km s ' (~ 
beam). Refer to Fig{T]for a guide to the symbols used for the compact objects. 
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Fig.A.3. SMA/APEX channel maps of the H2CO 322-221 line. Contours every 0.5 Jy beam ' km s ' 
beam). Refer to Fig{T]for a guide to the symbols used for the compact objects. 
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Fig. A.4. SMA/APEX channel maps of the CH3OH 422-3 12 line. Contours every 0.5 Jy beam ' km s ' (~ 3cr at edge of primary 
beam). Refer to Fig{T]for a guide to the symbols used for the compact objects. 
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Fig. A.5. SMA/APEX channel maps of the HC3N 24-23 line. Contours every 0.5 Jy beam ' km s ' (~ 3cr at edge of primary beam). 
Refer to Fig[T]for a guide to the symbols used for the compact objects. 
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Fig. A.6. SMA/APEX channel maps of the C-C3H2 6i6-5o5 and 6o6-5i5 lines (blended). Contours every 0.5 Jy beam"' km s"' (~ Bo- 
at edge of primary beam). Refer to Fig{T]for a guide to the symbols used for the compact objects. 
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Fig. A.7. SMA/APEX channel maps of the C-C3H2 5i4^23 line. Contours every 0.5 Jy beam"' km s"' (~ 3cr at edge of primary 
beam). Refer to Fig{T]for a guide to the symbols used for the compact objects. 
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Fig. A.8. SMA/APEX channel maps of the C-C3H2 524^13 line. Contours every 0.5 Jy beam"' km s"' (~ 3cr at edge of primary 
beam). Refer to Fig{T]for a guide to the symbols used for the compact objects. 
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Fig. A.9. SMA/APEX channel maps of the DCN 3-2 line. Contours every 0.5 Jy beam"' km s"' (~ 3cr at edge of primary beam). 
Refer to Fig[T]for a guide to the symbols used for the compact objects. 
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Fig. A.IO. SMA/APEX channel maps of the '^CN N = 2-1, J = 5/2-3/2 = 3-2, F = 4-3; = 3-2, F = 3-2; and Fi = 3-2, F 
= 2-1 lines. Contours every 0.5 Jy beam"' km s"' (~ 3cr at edge of primary beam). Refer to Fig{l]for a guide to the symbols used 
for the compact objects. 
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Fig. A.ll. SMA/APEX channel maps of the '^CN = 2-1, 7 = 3/2-1/2 Fi = 1-0, F = 0-1 line. Contours every 0.5 Jy beam"' km s"' 
(~ 3cr at edge of primary beam). Refer to Fig{T]for a guide to the symbols used for the compact objects. 
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Fig. A.12. SMA channel maps of the '^CN N= 2-1, J=5/2-3/2 F = 5/2-3/2, F = 7/2-5/2, and F = 3/2-1/2 lines. Contours every 
0.5 Jy beam ' km s ' (~ 3cr at edge of primary beam). No APEX short-spacing data were used for this line. Refer to Fig[l]for a 
guide to the symbols used for the compact objects. 
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Fig. A.13. SMA/APEX channel maps of the SiO 5^ Hne. Contours every 0.5 Jy beam ' km s ' (~ 3cr at edge of primary beam). 
Refer to Fig[T]for a guide to the symbols used for the compact objects. 
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